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Abstract 

We have been modeling the envelope of the 
secular light curve (Ferrn, 2010) [3 of comet 
lP/Halley. As the first step we have mod- 
eled the water production rate of the comet, 
for several active regions as a function of vi- 
sual albedo and infrared albedo, and then we 
have used the correlation between the water 
production rate and the reduced magnitude 
(Jorda, 2008) [5 for modeling the envelope 
of the secular light curve of the comet. We 
obtain probable orientations of the rotation 
axis (I,3>). These orientations are compared 
with several solutions by several authors. We 
have calculated the surface temperature for 
the orientation of the rotation axis with min- 
imal standar deviation. We have found that 
for near and far values of the perihelion the 
surface temperature as a function of the lat- 
itude is constant. 

1. Introduction 

Comet Halley has a large number of observed 
apparitions, and many visual observations 
have been published in the literature. Ferrm 
(2010) [3 has published brightness observa- 
tion of several comets, considering that the 
envelope of the data describes the behavior of 
the brightness in time. Several secular light 
curves show asimmetry. We have developed 
a model that predicts this asimmetry. The 
orientation of the rotation axis explains the 
asimmetry in the secular light curve (Rondon 
& Ferrn, 2010; Rondon, 2007)00. The im- 
portance of studying the secular light curve 
is that it gives a large amount of physical in- 
formations of comet. By modeling this curve 
we can explain the behavior of the brightness 



mathematically and we can predict the phys- 
ical parameters of the curve. 

2. Model Calculations 

The first step for modeling the light curve 
is to calculate the sublimation rate of water. 
The equation that describes the vaporization 
rate of a comet is the energy conservation 
equation, given by: 



F (l - A v )r^cos(0) = (1 - A ir )aT 4 + 

Z(T)L(T)+K^1 (1) 

Where A^ r is the infrared albedo, th Sun- 
comet distance, cos(0) is the projection factor 
(Cowan&A'Hearn, 1979) [T], a is the Steffan 
Bolztmann constant, T is the temperature, 
Z(T) is the sublimation function, L(T) is the 
latent heat function, K is the thermal con- 
ductivity constant, z is the layer depth. 
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2(tt)/cT 

where P(T) is the vapor pressure function, 
m is the molecular weight, k is the ideal gas 
constant. 

If we know the vapor pressure function and 
the latent heat function, we can solve for the 
energy conservation equation(l). 
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Ztotai = Z(i) = - (Z(i,b)cos(b)db) 

Z J-ir/2 

(3) 

We have modeled the water production 
rate for comet lP/Halley using the obser- 
vational data given by (Schleicher, 1998) [9 , 
and we have found that when considering an 



active region of 80 km we can predict the 
average of the observational data (Fig. la). 
However, the secular light curve considers the 
envelope of the data (Ferrin, 2005) [2 j. For 
this case we have taken an active region of 
180 km(Fig.lb). In both cases we have as- 
sumed an visual albedo A v = 0.009, an in- 
frared albedo Air = 0.5 and a thermal con- 
ductivity K=0. 

The correlation equation between the re- 
duce visual magnitude and the water pro- 
duction rate is known (Jorda, 2008) [5] and 
is given by: 

m = 125.051 - 4.077 \og(Z total ) (4) 




Obliquity(degree) 

Figure 2: Standard deviation for each of the 
orientation of the rotation axis for the comet 
lP/Halley. 



The secular light curve of comet Halley is 
assymetric. This can be explained through 
of the orientation of the rotation axis of 
comet (Rondon, 2007)(Rond6n & Ferrin, 
2010)00. 
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Figure 1: Water production rate vs log(r), a) 
with an active region of 80 km b) with an 
active region of 180 km. 



We have calculated the standar deviation 
for our model with the observational data, 
and found that the solution with minimal 
standar deviation is for a obliquity, I = 90 
and a pole orbital longitude, = 112, (Fig, 
[2J. Several authors obtained solutions for 
the orientation of the rotation axis (Sekan- 
ina, 1986) [TD]. 
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Figure 3: Secular light curve of the Halley 
comet with a orientation of the rotation axis 
of I = 90, = 112. In the upper envelope 
curve has been considered an active region of 
180 km and in the lower envelope curve has 
been considered an active region of 80 km. 

In Figure [2] we can see that a region of the 
solution of (Grun, 1986) [4 has a standard 
deviation of only 1 magnitude. In the Figure 
[3] we have modeling the secular light curve 
assuming an orientation of the rotation axis 



of 1= 90, $=112 for an active region of 80 
km(blue line), and 180 km, (red line). In this 
graph we can see that before perihelion the 
model fits the observational data, but after 
perihelion the fits is not very good. This 
behavior is due to the thermal conductivity 
effect in the nuclear surface that has not 
been considered in the model. In Figure 
4 we can see that the surface temperature 
as a function of the latitude is constant for 
values of the heliocentric distance near and 
very far away of perihelion. This behavior 
near of perihelion is that the coma of comet 
tends to stabilize the temperature, while for 
heliocentric distance very far away perihelion 
the sun radiation is small and therefore the 
whole cometary nucleus cools. 
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Figure 4: Surface temperature as a function 
of heliocentric distance, assuming an active 
region of 180 km and a orientation of the ro- 
tation axis of 1= 90, $ = 112. 



that the surface temperature as a function of 
the latitude is constant for values of the he- 
liocentric distance near and very far away of 
perihelion. 
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3. Summary and Conclu- 
sions 

We have developed a theoretical model capa- 
ble of reproducing the observational data of 
the secular light curve and water production 
rate of comet Halley, with minimal standar 
deviation a = 0.3 mag, for an orientation of 
the rotation axis 1 = 90, 0=112. The region 
of the solution shown by (Grun, 1986) [4] ex- 
plain the secular light curve with a standard 
deviation cr=0.6 5 mag. We also have found 



